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Abstract 

We develop a biorthogonal linear-scaling algorithm for a transcorrelated method based on the 
localized nature of transformed orbitals. The transcorrelated method, which employs a similarity- 
transformed Hamiltonian referred to as a transcorrelated Hamiltonian, enables highly accurate 
first-principles condensed-matter calculations in principle. Meanwhile, the transcorrelated Hamil- 
tonian numerically prevents us from applying it to large systems because the transcorrelated Hamil- 
tonian is a non-Hermitian operator and contains a 3-body electron-electron interaction term. Non- 
Hermiticity means that left and right wave functions of the total energy expectation value are 
different from each other. Namely, a biorthogonal form is required. Our new method allows us to 
handle the non-Hermitian operator and exhibits a linear-scaling behavior. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

The transcorrelated method , 0] , proposed by Boys and Handy, is re- 

garded as one of the approaches based on the Jastrow-Slater-type wave function, such as the 
Variational Monte Carlo 9|] method. The salient aspect of the transcorrelated method is its 
treatment of the Jastrow function. Instead of naively handling the Jastrow-Slater-type wave 
function, a similarity-transformed Hamiltonian with respect to the Jastrow function is used. 
Thanks to this particular Hamiltonian, referred to as a transcorrelated Hamiltonian, the 
many-body integral that stems from the function form of the Jastrow function is reduced to 
just a three-body one. Furthermore, we can obtain one-body energies that satisfy Koopmans 
theorem 7|. Considering that even within the single Slater determinant, the Jastrow-Slater- 
type wave function yields highly accurate ground-state energy, the transcorrelated method is 
thought to have great potential as a first-principles calculation tool. However, one must still 
tackle the non-Hermiticity of the transcorrelated Hamiltonian and the three-body electron- 
electron interaction term that is absent in the original Hamiltonian. The non-Hermiticity 
involves the biorthogonal treatment in which left and right eigenf unctions are different from 
each other. Furthermore, the three-body interaction term obviously limits the size of systems 
that we can investigate by the transcorrelated method, because it requires large computa- 
tional cost that scales as O {N^i)- Here, N^i denotes the number of electrons. 

The purpose of this paper is to propose an efficient O (A^ei) method of handling the above 
two issues with localized orbitals. In Sec. [Tll we introduce the transcorrelated method at 
first. And then we describe details of our approach in sec. IIIII In the first two subsections of 
Sec. we develop a biorthogonal self-consistent-field (SCF) equation for non-orthogonal 
orbitals. After that, we provide the important proof for handling the inverse of the overlap- 
ping matrix. Lastly, we discuss about the scaling behavior. 

II. TRANSCORRELATED METHOD 

At first, we briefiy introduce the biorthogonal transcorrelated method for the later dis- 
cussion. In the transcorrelated method, instead of the ordinary many-body Hamiltonian H, 
one addresses the transcorrelated Hamiltonian 

Htc = F~^HF. (1) 
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Here, F represents the Jastrow function. Because of the non-Hermiticity of the transcorre- 
lated Hamiltonian Htci the variational approach to the energy expectation value of HtC) 



E = I ^*HTc^d^''r 



(2) 



does not give the proper ground-state energy of the original Hamiltonian if, in principle. 
Instead, to obtain the proper ground state, the biorthogonal approach must be taken, in 
which left and right wave functions are different from each other, as 



(3) 



In the formulation given below, the transcorrelated method in this biorthogonal form is 
discussed with the single Slater determinant. Spin indices are omitted for simplicity. 
Suppose, for instance, the function form of the Jastrow function is given by 



exp 



1=1 



(4) 



then the transcorrelated Hamiltonian can be explicitly written as 



TC 



i=l 



1^2 



(5) 



with the notations 



i=l jj^i 
^ N^i N^i N^i 

g 

i=l j^i k^i,j 
1 



+ V^M (r-i,r,- 



- ViU {ri,rj) ■ ViU {ri,rj) 
Thus, the SCF equation and its Hamiltonian are written as 

HTC-SCF(f>Ri (x) = ^ eji(l)Rj (x) 



(6) 



(7) 



(8) 
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and 



HTC-SCF(t)Ri {x) 




(9) 



X \\(pRi{x)(l)Rj{y)\\ 




Here, {sij}, {(pLi}, {^ra}, and || - ■ ■ || denote constants originating with the Lagrange multi- 
phers, the left one-body wave functions, the right one-body wave functions, and the Slater 
determinant, respectively. It should be noted that the one-body wave functions {4>Li,4'Ri} 
are biorthogonal. 

III. LINEAR-SCALING METHOD IN BIORTHOGONAL FORM 

Here, we discuss a measure against the unfavorable scaling behavior on the number of 
electrons N^i. Utilizing the localized nature of transformed wave functions referred to as 
localized orbitals, the scaling behavior can be improved from O (N^i) to 0{Nf.i). Yet the 
localized orbital is not a kind of orthonormal wave function in general; thus, Eq.® should be 
replaced with one without an orthogonality restriction. This idea is regarded as an extension 



of the work of Mauri et al. Il2l . 




for linear scaling calculations. 



A. Energy expression without orthogonality restriction 



Let us consider, for instance, the two-body interaction term 



Tr [vee (x, x') p {x, x'; y, y')] 



(10) 



where 




(11) 



X F {x, x, xs,-- - , XnJ ^* {y, y', ajg, • ■ ■ 



XnJF 
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with ^ being the exact many-body wave function of the ordinary many-body Hamilto- 
nian H. The two-body reduced density matrix Eq. ffTTl) has a different definition compared 
with the ordinary one H due to the existeuee of the Jast.ow funct.on F The^efo.. 
here and hereafter, we call it the similarity-transformed two-body reduced matrix. As dis- 
cussed in Appendix, since left and right eigenvectors of the transcorrelated Hamiltonian Htci 
{^F,F~^'^}, in Eg. (ITT]) are approximately treated as Slater determinants in the transcor- 
related method, the similarity-transformed one-body reduced density matrix, 

J (12) 
X F"^^ {x, X2,-- - , xnJ ^* {y, X2,-- - , XnJ F, 

satisfies the idempotency condition as in the case of Hartree-Fock method, and the cumulant 
expansion is also applicable to the similarity-transformed two-body reduced density matrix 
Eg. (ITT!) as follows : 

P {x, x'; y, y') = {x; y) p {x'] y') - {x'] y) p [x] y') . (13) 

Also, according to the proof in Appendix, the similarity-transformed one-body reduced 
density matrix is expressed in terms of its left and right eigenvectors, {ipLi} and {ipRi}, as 

p{x-y) = Y,i'Rd^)ru{y)- (14) 

i=l 

Consider the non-biorthonormal orbitals {(pu} and {(pm} defined as 



N, 



el 



hi {x) = ^ UJijipLj (x) 



{x) = ^ uJijipRj (x) . 



Conversely, from Eg. (fT5|) . one obtains 



i^Li (X) = J2^ij^^Lj{x) 



(16) 



i'Rii.x) = ^UJ-^^(t)Rj{x) . 
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Now, the overlapping matrix between {(pLi} and {4>Ri} is 

Sjk = j dxcplj (x) (f)Rk (x) 



^ uj*^ujki' j dxi^li (x) iJni' (x) 
,i'=i ^ 



i=l 



ik ■ 



Therefore, the inverse of the overlapping matrix is 

i=l 

Note that tu^^* means (tu"^)*^. Substituting Eqs. f[T6l) and (|T8l) into Eq. (|T4l) yields 



p {x; y) = Yl '^^^ ^y") 



1=1 



i=l lj=l 



Nei 



'Nel 
.k=l 



Nel Nel ( Nel \ 

fc=l j=l \i=l J 



Nel Nel 

= E E (^) 02. (y) ■ 

k=i j=i 

As a consequence, the two-body interaction term becomes 
Tr [vee {x, x') p (x, x' ] y, y')] 
= J dxdx'dydy'5 {x — y) 5 {x' — y') v^.^. {x, x') 
'I 1 

X 2^ p y'^ ~ 2^ y^ ^ y'^ 

= J dxdx'dydy'5 {x — y) 5 {x' — y') v^e (x, x') 

Nel 

X2 E S7^S^,^[<PRdx)<Pl,{y)<pRk{x')<Pl,{y') 

i,j,k,l=l 

- ct>m {x') 02, {y) cl>Rk (x) <t>l^ iy') ] 



. Nel I- 

2 E Sr^S,n dxdx'<Pl^{x)<Pl,{x')v,,{x,x') 



i,j,k,l=l 
X \\(l)Ri{x)^Rk{x')\\ 



In a similar fashion, since the three-body reduced density matrix can be decomposed as 



p (33, x', x"; y, y', y") 



p{x;y) pix;y') p{x;y") 
p{x';y) p{x'-y') p{x';y") 



p{x";y) p{x";y') p{x";y"] 



the total energy without the orthogonality restriction is expressed by 

1, 



i=l 



+ ^ Vee (Xi, Xj) p (Xi, Xj] y^, y j) 

+ "^eee {Xi, Xj, Xk) p {xi, Xj, Xk] y^, yj, yj}j | 

i,j,k=l 

= Sij" / dx<Pl^ {x) 
i.j=i 

+ 2 5Z Sr-^S^l dxdx'<Plj{x)<Pli{x')v,,{x,x') 



--V'^ +Vext {x) 



^Ri [X) 



i,j,k,l=l 

X \\(f)ni (x) (pRk {x 

N. 



(21) 



(22) 



-1 JVel „ 

+ 6 E S^^SMS^n J dxdx'dx"4>l^ {x) 02, {x') 

i,j,k,l,m,n=l 

X 02„ {x") Veee {x, x' , x") ix) (j)Rk {x') (pRm {x")\\ . 



B. Biorthogonal SCF equation without orthogonality constraint 

The biorthogonal SCF equation is derived by the variational approach to the total energy 
E in Eg. (12^ . namely, dE/dcj)}^^ (w) = and dE/dcj)*^^ (w) = 0. For instance, the functional 
derivative of the two-body term E2 in the total energy E with respect to i'w) is 
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77^ = o P>^* '^f.,\ ^ki^ / d^dyi^lj (^) 'i'Li iy) ^ee (a;, y) Urn {x) <i)Rk (y) \\ 

^ i,j,k,i=i ^ > 

+ 2 E ^^^' d6* '(w) / '^^^ "^^^ ll*^^^ '^^^ ^^^11 

i,j,kl=\ ^ ' 

+ 2 XI ^ij^^k^ f d^dy^jp^ - (f^Ll (y) Vee {x, y) {x) (j)Rk (t/) || 

i,j,k,l=l 

+ 2 X '^ij^'^fe/ / dxdy(t)lj {x) Sip6 {y - w) v^e {x, y) (x) (t)Rk (y) \\ 

ij,k,l=l 

^ ,\ / ^^^Lj (^) 1 o ^kl / dycf^li (y) [Vee (x, y) + V^e (y, x)] 



X \\(f>Riix)(f)Rkiy 



^Y^iP^\\Y ^kl I dy4>*Ll (y) bee (w, y) + "Uee (y, 0iife {y 

i=l I fe,i=l 

^ Yl / ^^'^ij H^C-SCF2<pRi {X) + X] S¥^TC-SCF2<pRi M . 

(23) 

Here, H^fj gQp^ is the operator defined as 

HTC-SCF2<t>Ri (^) = 2 X "^^'^ / ^^-^ ^^^^ ^ "^^"^ ^''^ ^^-^ '^'^^ ■ ^^^'^ 

fc,i=l 

In exactly the same way, one readily obtains 

dE dS~^ f 

d^l^ ^ Y d^J^ J HTc-scF4>m {x) + Si-p'H^c-scF<Pm M , 

(25) 

where 

HTC-SCF4>Ri [x) 



-^V^ + Vext (X) 



(t>Ri + 2 XI ^Jk j dy(t>*Lk (y) bee (x, y) + V^e (?/, x)] {x) (t)Rj (?/) 

j,k=l 

+ g X ^Jk^hl j dydZ(j)lk {y) (j)*L^ {Z) [Veee {x, y, Z) + Veee (?/, Z, x) + Veee {z, X, y) ] 



j,k,l,m=l 

X \\(l)Ri{x)(l)Rj{y)(l)Ri{z)\\. 



(26) 
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The functional derivative of the inverse matrix 5*"^ with respect to 0^^^ (it>) in Eg. (1250 is 
not straightforwardly computable so we therefore compute it by utilizing the relationship 

SS^^ = I as, 



o-i c-i 
d(t>*Lp M 



ij 



(27) 



Substituting the above-mentioned result into Eq. fl25|) yields 
dE 



; 1 A — 1 

(28) 



^ ' i,j,l=l j=l 



Also, operating J2p=i ^pk from its left side, one has 



5Z '^Pfc a^* X = - Yl{Yl^l^j^ / '^^^Lj {X) HTC-SCF<pRk (X) \ (pRl (W) + H^^_scF<pRk [w) 

P=i ^"Plp^^) J J 

E^fflfe^ra (w) + H^c-scF^'Rk i'w) , 



p=l -rnpy I 



1=1 



(29) 



where 



= E ^ik / dy<Pl, {y) H^c-scF^m {v) ■ (30) 
k=i 

Because dE/dcjfj^p {w) = 0, the SCF equation is finally derived as 

Hrc-scF^pRi (x) = E ^RJi^RJ (^) • (31) 
i=i 

Similarly, H:^^_gfjp, the SCF equation, and eiji corresponding to {4>Li} are obtained as 
Eqs.gE^, dSH), and ([3D]), respectively. 

C. Handling S'^ 

The inve.e of the overiapping n,atrix. S-. .e,ui.e. O WJ) computation. It should 

therefore be replaced with the polynomial function of 5*1111. |12| , namely, 

Ns 

S-^ = Y,iI - ST (32) 

n=0 
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with an odd number Ns- As discussed later, it is straightforward to show that the com- 
putational cost for performing Eq. fl3T|) is reduced to O {Nei) by utilizing localized orbitals 
accompanied with Eg. (15^ . 

Here, we discuss whether the replacement ( 1^2]) causes any errors in the total energy. To 
do so, we decompose the overlapping matrix S with its eigenvalue matrix A, left eigenvector 
matrix V, and right eigenvector matrix U as 



S = UAV. 



(33) 



It is easily shown, by using 6 = "^^=0 (-^ ~ — ^ A)'^^^^ — I , that the eigenvalue 

matrix of Q — S^^ is reduced to — A^^ (/ + S) as 

v{Q- s-^) u = -A-^ + e 

= -A-^ [/ - AG] 

= -A-^ [/+(/- A) - 9] 



-A-i 



-A-i 



Ns+l 

1+ j2 i^-^y 

n=l 



e 



(34) 



/+(/-A) 



Ns+l 



= -A-' (J + S). 

The non-negative definiteness of S makes A~^ a non-negative definite matrix. Also, J + H = 
(/ — A)^'^'*'^ is a non- negative definite matrix for odd number Ns- The matrix Q — is 
therefore a non-positive definite matrix because both A~^ and J + H are non- negative definite 
matrices. Using A, U, and V, one obtains 

QijQkl = j UipQppVpj 1 I UkqQqqVgl j 

^ (35) 



.9=1 



UipUkqQppQgqVqlVpk. 

p,q=l 



This is a spectrum expansion of the N^i x N^^ matrix '^Q which has an element QijQki in the 
{i, k)-th row and the (j, /)-th column. Thus, QppQqq is the (p, g)-th eigenvalue because of 

QijQklUjpUlq = QijUjp QklUlq 



j = l 1=1 
0pp ® 59 ^ip Ukq ■ 



(36) 
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If Nil X Nil matrix ^S* Ms defined as 



el el 



(37) 



Q — S is proved to be a non-positive definite matrix as follows: 



6 6 - A'^A'^ 

~ ~^pp ^qq ~ ^pp^qq^pp^qq) 



-KiKq 



l-(-A0)pp(-A0) 



/qq 



(38) 



-A-'^A-^ (1 

pp qq \ 



■^PP^qqJ 



< 



within the assumption of < App < 2. By the same token, if N^^ x N^^ matrices and 



■^5* ^ are defined as 



and 



= {QijQkiQ 



el el 



i^ij ^kl ^mn) 



ei ei 



(39) 
(40) 



respectively, Q — S is also proved to be a non-positive definite matrix as 

6 6 6 - A~^A~^A~^ 

'^pp '^qq^rr ^pp ^qq ^^rr 

Aj,y. (1 AppAqqAjj.0pp0qg0j.fl) 

= -A'^A'^A'^ (1 

pp gg rr \ 



(41) 



■^pp' — 'qq' — 'rr ) 



< 0. 



As a result, according to the proof by Mauri et al. UJ, ll2| , it can be proved that even when 
is replaced with Q, one obtains the same total energy E by restricting App in the range 
[0, 2]. The restriction of App in this range is easily achieved by adjusting the trace of S. 



D. Scaling behavior 



Now, we discuss the scaling behavior of the above-mentioned scheme. For example, 
consider the exchange term in the total energy replaced with Q, 



2 XI <5iiQfc« / dxdy(f)lj{x)(j)li{y)vee{x,y) 

/, 1 1 



ij,k^l=l 

X (pRi {y) <pRk (x) 



(42) 
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The term Qij means that the index j runs over the adjacent indices to the index i. Similarly, 
Qki and in) 4>Ri iv) restricts the indices k and I on the adjacent indices to the index i. 
Ultimately, the only index i runs over all indices. Thus, this exchange term shows an O {N^i) 
behavior accompanied with localized orbitals. 

Next, consider the three-body interaction terms. Classifying by the number of inter- 
changed pairs, there are three cases, namely, for example. 



(43) 



(44) 



(45) 



Et"" = ^ QijQklQmn J dxdx'dx"(t>l^{xWu{x') 

i,j,k,l,m,n=l 

X (pLni^")'"eee{x, x' , x")(j)Ri{x)(j)Rk{x')(j)Rm{x") , 

1 r 

=6 ^ QijQklQmn I dxdx'dx"ct>l^{x)ct>li{x') 

i,j,k,l,m,n=l 

X <f>Lni^")Veee{x, x' , x")(l)Ri{x)(f)Rm{x')(f)Rk{x") , 

1 r 

=6 ^ QijQklQmn J dxdx'dx"(t>l^{x)(t>li{x') 

i,j,k,l,m,n=l 

X (pLni^")Veee{x, x' , x")(f)Rm{x)(f)Ri{x')(f)Rk{x") . 

In the non-interchanged case, E^°'^, one obtains 

-, N^i ^ 
Es""" = g QijQklQmn J dxdx'dx"<Pl^{x)(Pli{x') 

i,j,k,l,m,n=l 

X (l)*Ln{x")Veeeix, x' , x")(f)Ri{x)(f)Rk{x')(f)R^{x") 

^\J^^\Y1 Qij'pLjix)(l>Ri{x) } J dx' <Y Qki4>li{x')(t>Rk{x') ^ (46) 

X / rfa;Vee(a;,Cc',a3") < Yj Qmn<P*Lnix")(pRmix") 

ym,n=l 

= l/da=p(.;a=)/d.V(.';.')/<i.V„(.,. ',.>(.»;."). 

The cost which we pay for this is just an 0{Nei) to calculate p{x;x). Therefore, it is 
computable with an O^N^i) cost. Similarly, the one-pair interchanged case, E^"^^, can be 
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rewritten as 



=1 E QijQklQmn j dxdx'dx"(t>l^{x)(t>li{x') 

i,j,k,l,m,n=l 

X (pLni^")'"eee{x, x' , x")(f)Ri{x)(f)Rm{x')(f)Rk{x") 

^ 6 / I XI ^^-J'^^J'^^)^^^^^) f XI QklQmn j dx'dx"(j)li{x') 

(.i,j=l J k,l,m,n=l 

X 02„(a;")^eee(a;, x' , x")(i)Rm{x')(t)Rk{x") 

dxpi.^.) Y^ Q^^jd^d^'^.i^) 



(47) 



1 

6, 

fe,/,m,n=l 

X (i)l^{x")veee{x, x! , a;")0R^(cc')0Rfc(cc")- 

As in the case of the exchange term, only one index among k,l,m, and n runs over all indices. 
Thus, it also shows a linear-scaling behavior. Note that p{x; x) can be computed separately 
and beforehand. In the case of the two-pair interchanged case, Ej'^", the terms Qij, Q^u 
Qmn, (l)ljix)(pRrn{x), (j)li{x') (j) Ri{x' ) , and (j)\,^^{x")(j)Rk{x") restrict 5 indices. Thus, only one 
index among 6 indices runs over all indices. As a result, the three-body interaction terms 
also show a linear-scaling behavior. Likewise, the other terms which are not mentioned here 
in the total energy and the SCF equation are computable in O {N^i) computational costs. 



IV. CONCLUSION 

In summary, we have proposed an efficient approach for the transcorrelated method in the 
biorthogonal form. In this approach, the original transcorrelated Hamiltonian is replaced 
with the one free of the orthogonality constraint in order to handle the non-orthogonal 
localized orbitals. Because of being free of the orthogonality restriction, the unhandy inverse 
matrix S"^ appears in the formulation. To achieve a linear-scaling behavior of the whole 
framework, we replaced S"^ with the polynomial function of S referred to as Q, and we 
gave the proof that even when is replaced with Q, the same ground state is obtainable. 
Lastly, we showed that our approach has an O {N^i) scaling behavior. 

For now, in this paper, we proposed and discussed the theoretical framework only. To 
prove the effectivity of our approach, a number of numerical examples are required indeed. 
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Currently we are under investigation in this direction. 
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APPENDIX A: CUMULANT EXPANSION OF THE SIMILARITY- 
TRANSFORMED DENSITY MATRIX 

Consider the exact many-body wave function which satisfies the shrodinger equation 

= E^. (Al) 

Then the exact left and right eigenvectors of the transcorrelated Hamiltonian 
Htc{= F~^HF) are written with this * as 

*L = *F, (A2) 

= F-^*. (A3) 

These two eigenvectors are approximately treated as Slater determinants in the transcorre- 
lated method, namely, 

= \\'4^R1^R2 ■ ■ ■ i^RN^i II • (A5) 

We are not approximating the exact many-body wave function ^ solely. We have a high 
degree of freedom for {^Lh'^Ri} as in the case of the Hartree-Fock method, because Slater 
determinants are invariant under any unitary transformations. So we can select biorthogonal 
sets as {'ipLi, ^ffi}- From the things mentioned above, the similarity-transformed many-body 
density matrix is derived as 
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p{xi,x2,--- ,XN,i;yi,y2, - ■ ■ ,yNj 
= ^r{xi,X2,--- ,a;jv,J^l(z/i,2/2,--- ^VnJ 



Nell 



Nell 



^m(a^2) ^R2{x2) 

p{xi;yi) p{xi;y2) 
p{x2;yi) p{x2;y2) 



^Ll(2/2) V'L2(Z/2) 



i^LNAyi) 

'^LN,Sy2) 



^LiiywJ ^L2{yNj ■■■ i^LN^iiywJ 



pi^'i-^yNj 



pixN,^;yi) p{xN,{,y2) ■■■ p{xN,,;yNj 



Here 



p{x;y) = '^ijRi{x)ijli{y). 



(A6) 



(A7) 



i=l 



This is exactly the same definition with that of the reference 15| except the definition of 
the one-body reduced density matrix. This difference between definitions of the one-body 
one leaves no effect on the applicability of the cumulant expansion when the one-body one 
is still idempotent. In other words, if the one-body one satisfies 



dyp{x; y)p{y; x') = p{x; x'), 



(A8) 



the cumulant expansion is also available to the similarity-transformed ones. It is readily 
shown that the biorthogonality of {ipLi, i^Ri} proves an idempotence of the one-body reduced 
density matrix. 



a. Proof of the availability of the cumulant expansion to the similarity-transformed reduced 
density matrix 

At first, we define a similarity-transformed M-body reduced density matrix as 
p{xi, X2,--- , xm; 2/1, ?/2, ■ ■ ■ , 2/Af) 



[ j^j r^Rixi,X2,--- ,XNj'^l{yi,y2,--- ,yNj H dxkdyj,5{xk-y 

^ ^ k=M+l 



(A9) 
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From this definition, a A'^erbody one is 
p{xi, X2, ■ ■ ■ , a^Tv^p 2/1, 2/2, ■ ■ ■ , VnJ 

p{xi\y{) p{xi;y2) ■■■ p{xi;y^J 

p{x2]yi) p{x2]y2) ■■■ p{x2]yNj 



Nell 



(AlO) 



pixN,,;yi) p(a;jv,n2/2) ■•• p(a=JV,,; ^tv.J 

as shown in Eg. ( ]A6[) . Also it is easily checked that a (M — l)-body one can be obtained 
from a M-body one as 



p{xi,X2,--- , Xm^i; 2/1, 2/2, ■ ■ ■ , 2/A/-1) 

\m-i) 



\M ) 



M 



Nel-M + 



dxMdyM^XM - yM)p{xu 2:2, ■ ■ ■ , Xm] 2/1, 2/2^ ■ ■ , 2/m) 
Y j dxudy^SixM - yj^)p{xi, a;2, • ■ ■ , a^A/; 2/1, 2/2, - ■■ , 2/m)- 



(All) 



When the M-body one is given as 

p{xi,X2,--- ,a;A/; 2/1,2/2^ ■ ■ ,2/ A/) 

p(a3i;2/i) p{xi;y2) ■■■ p{xi;yM) 
p{x2;yi) p{x2;y2) ■■■ p{x2;yM) 



1 

M! 



p{xM]yi) p{xM;y2) ■■■ pixM;yM) 



(A12) 



the (M — l)-body one is 



p{xi,x2, ■ ■ ■ ,XM-i;yi,y2, ■ ■ ■ .Vm-i) 
M 



N,i-M + 



- j dxMdyj^.j5{xM - yM)p{xi^x2, - ■ ■ ,XM;yi,y2, ■ ■ ■ ,yM) 



{M-iy. 



dxMdyM^ixM - yu) 



P{xi\y^) p{xi;y^) ■■■ p{xi;yj^^) 
p{x2;yi) p{x2;y2) ■■■ p{x2;yM) 

p{xM;yi) p{xM;y2) ■■■ p{xM;yM) 
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{N,i-M+1)-^ 

{M-iy. 



dx 



M 



{M-iy. 



M-l 



p{x\\Vi) P{xi\y2) ■■■ p{xi;xm) 

p{xM;y-i) p{xM;y2) ■■■ p{xm;xm) 

p{xi;y{) pixi;y2) ■■■ p{xi;yM-i) 

, , , pix2\y^ p{x2;y2) ■■■ P[X2;yM-ij 
dxM\ p{xM;yM) 

p{xM-i;yi) p(xm-i;2/2) ••• pixM-i;yM-i) 

P{xi;yi) ■■■ p{xi;y^_^) p{xi;yk+i) ■■■ p{xi;xm) 
p{x2]yi) ■■■ p{x2;yk-i) p{x2]yk+i) ■■■ p{x2;xm) 



fe=i 



(N,, -Ad + iy 

(M - 1)! 



M-l 



+ E (-1) 



M+k 



fe=l 



p(xm-i;2/i) ■•• p{xM-i;yk-i) p{xM-i;yk+i) ■■■ p{xm-i\xm) 

p{xi;yi) p{xi;y2) ■■■ p{xr,yM-i 

^ p{x2;yi) p{x2\y2) ■■■ p{x2;yM-i) 
dxMP{xM;yM) 

p{xM-i;yi) p(xm-i;2/2) ••• p{xM-i;yM-i) 

P{xi;yi) ■■■ p{xi]y^_^) p(xi]y^.^^) ■■■ J dxMp{xi;xM)p{xM;yk) 

p{x2;yi) ■■■ p{x2]yk-i) p{x2;yk+i) ■■■ J dxMp{x2\XM)p{xM;yk) 



p{xM-i;yi) ■■■ p{xM-i;yk-i) p{xM-i;yk+i) ■■■ J dxMp{xM-i;xM)p{xM\yk) 

P{xi;yi) p{xi;y2) ■■■ p{xi;yM-i 

(iVe; -M + 1)"^ f p{x2;y-i) p{x2;y2) ■■■ p{x2\yM-i) 



(M-l)! 



Nel 



M-l 



(^_l-^M+k^_^-^M-k-i 



p{xM-i;yi) p{xM-i\y2) ■■■ p{xM-i;yM-i) 

P{xi;yi) p{xi;y2) ■■■ p{xi;yj^_^) 

p{x2]yi) p{x2;y2) ■■■ p{x2;yM-i) 



fe=i 



p{xM-i;yi) p{xM-i\y2) ■■■ p{xM-i;yM-i) 
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(M-iy. 



(M - 1)! 



p{xM-i\yi) p{xM-i]y2 
p{xi;yi) p{xi;y2) ■■■ p{xi;y^.j_i) 
p{x2]yi) p{x2;y2) ■■■ p{x2;yM-i) 



p(a;i;yA/-i) 

P{x2]yM-l) 

p{xM^i]yM-i) 



p{xM-i;yi) p{xM-i;y2) ■■■ p{xM~i;yM~i) 

As a consequence, by the mathematical induction, if the A^^erbody density matrix is given 
as in Eg. flA6P and the one-body reduced density matrix is idempotent, then the M-body 
reduced density matrix for an arbitrary M G [1, A'e;] is expressed as 



p{xi;yi) p{xi]y^) ■■■ p{xi;yM) 
p{x2;yi) p{x2]y2) ■■■ p{x2;yM) 



1 
M! 



p{xM]yi) p{xM;y2) ■■■ p{xM;y 



Mj 



(A13) 
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